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Abstract ARTICLEINFO 

Digital image processing plays a pivotal role in various fields, Artie history: 

from medical imaging to surveillance systems. However, the Received 13 Jun 2023 

acquired images are often susceptible to various types of Revised form 15 Jul 2023 

noise, such in the form of salt and pepper noise, which can Accepted 25 Aug 2023 

severely degrade image quality and hinder subsequent 
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In this study, we introduce a fuzzy impulse noise removal es wi Fuzzy logic, Noise 
etection. 


algorithm as a potential solution. The efficiency of the 
suggested algorithms is assessed by comparing their 
performance to various existing noise removal methods. 
Objective measurements, including peak signal-to-noise ratio 
and mean square error, are used to evaluate the results. The 
findings demonstrate that the proposed algorithms deliver 
excellent outcomes in noise reduction and image detail 
preservation across a broad range of noise densities. 
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1. INTRODUCTION 

Noise can be described as any undesirable signal that contaminates an information-carrying signal It exists to 
varying degrees in almost all environments. Images often suffer from noise, which can emerge during the 
image acquisition process, transmission, or reproduction. The removal of noise from images is a crucial task 
in image processing and represents a significant obstacle to effective image analysis [1-6]. 


Among the different noise types, impulse noise holds significant prevalence and impact in digital images. 
This form of noise consists of brief, sporadic instances of "on/off" noise pulses that can introduce 
disturbances into images during the acquisition stage due to sensor-related issues (such as switching 
fluctuations or sensor temperature variations), or during transmission caused by channel irregularities (such 
as interference or atmospheric disruptions)[7-13]. Additionally, impulse noise can emanate from flawed 
hardware memory positions or synchronization instances discrepancies (similar to errors encountered during 
analog-to-digital conversion) that occur during the image processing phase [14,15]. 


and also Salt and Pepper noise, a type of image distortion, adds random black and white pixels to an image, 
resembling grains of salt and pepper. It can degrade image quality but is commonly encountered in real- 
world scenarios, requiring effective noise reduction techniques for restoration.[16,17] 
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Zadeh introduced fuzzy logic in 1965 as a mathematical tool designed to handle uncertainty. It presents a 
valuable approach for soft computing by incorporating the notion of computing with words and offering a 
method to address imprecision and information granularity. Within the framework of fuzzy theory, linguistic 
constructs like "many," "low," "large," "dark," "bright," and others can be effectively represented [18,19] 


2. Proposed work 

The proposed work comprises two primary stages. Firstly, it involves fuzzy noise detection to identify noisy 
pixels using suitable fuzzy sets and rules. Each model of impulse noise is treated differently in this stage 
based on the outcome of the noise identification process. Secondly, the work includes fuzzy noise filtering to 
recover the pixel values classified as noisy during the identifying noise stage. This filtering is achieved using 
appropriate fuzzy sets and rules as well. 


2.1. Impulse Noise Detection 

In numerous scenarios, applying impulse noise filtering without first classifying the noisy and noise-free 
pixels can lead to image distortion and blurring of edges. As a result, impulse noise detection becomes a 
crucial preliminary step before proceeding with noise filtering. 


Impulse noise is commonly caused by bit errors during data transfer and is characterized by its random and 
sparse nature. It stays autonomous and uncorrelated with the pixel elements of the image, resulting in some 
image pixels being noisy while the rest remain noise-free. The impulse noise model can be succinctly 
described as follows [20-24]: 


n(x, y), with probability P 


I(x,y) = aes (1) 
O(x,y), with probability 1 — P 


In this context, let I(x, y) represent the pixel value of the noisy image, O(x, y) denote the uncorrupted 
(original) image pixel value, and n(x, y) depict the pixel with impulsive noise at the coordinates (x, y). 


In this study, the impulse noise considered is Salt and Pepper noise (SPN). In Salt and Pepper noise, the 
noisy pixels can take on either the maximum intensity value L_max (gray level 255) or the minimum 
intensity value L_min (gray level 0). Consequently, the noise appears as white and black spots on the 
images. As a result, the value of n(x, y) in Equation (1) will be either 255 or 0. If (P) is the total noise 
density, then the noise density for salt noise and pepper noise will each be (P/2). In some cases, Salt noise 
and pepper noise could exhibit dissimilar levels of noise densities, denoted as Pl and P2, respectively, and 
the cumulative noise density will be given by: P= P1 + P2. 


When an image is affected by Salt and Pepper noise (SPN), a noisy pixel assumes one of two extreme 
values, either 0 or 255. Consequently, a pixel having a value of 0 or 255 is classified as a suspected pixel. 
However, if a suspected pixel closely resembles devoid of noise neighboring pixels, it can be considered as 
devoid of noise too. Two pixels are deemed similar when their absolute difference in grey value is small but 
not zero. This is because, in cases where both pixels are noisy, their absolute difference can be equal to zero. 


The term "small" to describe the absolute difference can be represented by a fuzzy set called "Small 
Absolute Difference," with its membership function p_ small illustrated in figure (1) 


Usmall 


“Small Absolute 
Difference” 


0 = 5 Difference 


Figure (1): Fuzzy set "Minor Absolute Deviation'' membership function 
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The membership function p for "small"is defined by two predefined parameters, namely a and b, with values 
a=10 and b=22 as established in [25-30]. 


The essential steps for detecting Salt and Pepper Noise (SPN) can be outlined as follows: 


Consider a neighborhood of size (2K+1)x(2K+1) around a central pixel I(x, y) at position (x, y) in the image 
I, where K starts with 1. If the value of the central pixel within this window does not match either of the Salt 
and Pepper noise values (i.e., I(x, y) # 0 or 255), it is identified as a pixel without noise. In any other case, 
the following steps will be executed: 


1. To determine the count of noise-free pixels within the observed window, we utilize the binary matrix 
M_SP as given in Eq. (3.2) along with the following equation: 


K K 
GZ, = >, py Mgp(x + s,y + t) 
Se sen (2) 

If ( GL, < 1), then the 

The window size will expand by increasing the value of K. This process is reiterated until the condition is 
met.(G¥, = 1) is met. 

2. Compute the absolute variances involving the central pixel and its neighboring counterparts within the 

specified window. as described below: 


Dyy =|I(x + s,y + t) —1(x,y)|, with I(x + s,y + t) # I(x, y) M3) 
where: s,t € {-K,..., +K} 

3. Find the maximum value in D,,, as: 

Myy = max(D,y) .- (4) 


4. The membership function Ls mq is used to determine whether the value of m,, is small as given in the 
following equation: 


i 1<my<a 
b-m 
Renan) = Tera as Myy <b w (5) 
0, Myy > b 


5. A fuzzy set named "Noise-Free" is employed to ascertain whether the current pixel I(x, y) can be 
classified as noise-free. The membership function for this set is derived as follows: 


Usmatt (Mxy) if I(x, y) = (0 or 255) 

Enoisefree U(x, y)) = 22 (0) 
1 otherwise. 

2.2. Filtering Impulse Noise 
The filtering procedure is applied exclusively to pixels with a membership degree less than one in the Noise- 
Free fuzzy [22][23] set. (LegUnowepacet (x,y)) <1 is considered as a noisy pixel). The filtering window 
size (2L+1)x(2L+1) is chosen adaptively according to the output of noise detection stage 
(1.€., Unoisefree) a8 follows: 


1. Calculate the count of noise-free pixels within that window by: 


Gy = a) Dia-11 (x +s,y + t) with leatsepresl IO +S,y + t)) =1 ...(7) 
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Di: AE (Gh, < 1), Subsequently, the window size will be enlarged by increasing the value of L. This process 


is reiterated until the specified condition is met (Gy > 1) is met. 


In cases where neighboring pixels display similarity or pertain to a notably uniform area, median-based 
algorithms prove successful in eliminating impulse noise. Nonetheless, in scenarios where pixels are part of 
regions encompassing edges or intricate image intricacies, the integration of information from other pixels, 
as represented by their corresponding weights, becomes imperative. Consequently, For every examined 
filtering window with dimensions (2L+1)x(2L+1), a fuzzy set denoted as 'Similar' is formulated to ascertain 
the level of similarity for each pixel within the window. This set is then utilized to allocate the suitable 
weight during the noise detection phase. 


Pixels demonstrating akin intensities within the observed filtering window will display a heightened 
membership degree within the 'Similar' fuzzy set. Likewise, edge pixels, characterized by comparable 
intensities along the edges, will also exhibit elevated membership degrees. Conversely, pixels affected by 
noise, deviating in intensities from their adjacent pixels, will manifest lower membership degrees. The 
‘Similar’ fuzzy set is depicted through a Gaussian-shaped membership function denoted as p similar, 
visually represented in figure (2). 


Hsimilar 


“Similar” 


Pixel intensity 


oO c 255 


Figure (2) Fuzzy set "Similar" characterized by a membership function resembling p_similar. 


The membership function for the fuzzy set "Similar" is defined through two parameters, namely, c and o. 
The parameter c designates the center point of the membership function, where the maximum value is 
attained. Meanwhile, o is associated with the distribution width of the membership function. As a result, 
these parameters are dynamically determined based on the homogeneity within the neighborhood of size 
(2L4+1)x(2L+1), as outlined below: 


c= mean(I(x + s,y + t)), with Posed OS ytOH=1 ... (8) 
—Lss,t< 
. 9) 


o = max (mean(|I(x +s,y +t) —cl), 0.01) 
—Lss,tsL 


The behavior of the membership function L,jmjiqgr can be described depending on Eq. (8) and (9) as follows: 


1. The extent of the membership function's width will be relatively greater in cases where the observed 
window belongs to a non-homogeneous region, such as edges or image textures. Conversely, the width 
of the membership function will decrease as the homogeneity level increases. In rare instances where the 
noise-free pixels within the observed window have identical values, the value of o will be zero, 
potentially leading to an infinite result. To avoid this situation, the Max operation in Eq. (9) is used with 
a small value (e.g., 0.01) to prevent encountering infinity in Eq. (10). 


2. The membership function reaches its maximum level at the average value of noise-free pixels within the 
observed window. 
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Consequently, by utilizing the membership function p similar, the effectiveness of the filtering process is 
enhanced. This leads to the removal of noisy pixels while endeavoring to preserve image details to the 
greatest extent possible. The expression for the membership function similar is as follows: 


Lae 


Usimitar Uy) = o( - ... (10) 


The final fuzzy weight w; for each pixel p, in the observed window of size (2L + 1) x (2L +1) is 
determined by the following fuzzy rule: 


Fuzzy Rule Establishing the degree of fuzzy weighting forp,: 
IF (p;, is noise-free) AND (py, is similar) 
THEN (w,, is high) 
This rule can be executed through the intersection operation involving two fuzzy sets. Consequently, the 
validity of the aforementioned rule is achieved by: 
Wz = Min{Unoisefree (Px), similar (Px)} . (1) 
In this context, the index k ranges from | to (2K+1)‘%2, enabling the selection of a window element. 


Ultimately, the outcome of the fuzzy filtering process for a specific pixel I(x,y) within the analyzed window 
of dimensions (2K+1)x(2K+1) is denoted as F(x,y). The calculation is performed as follows[6]: 
ap Vo ple ee ew Sree 
ee K Lee KWX+s,y +t) 
Here, /(x + s,y +t) denotes the pixels within the examined window surrounding the central pixel, while 
(x + s,y + t) signifies the associated weight assigned to each pixel within that window. 


+ w(x, y). 1 (a b2) 


3. Experimental Results 


In this section, The effectiveness of the suggested grayscale algorithm is assessed through a comparison with 
the following noise reduction methods: MF[11], HFF[12], AMF[13], FSM[14], and AWTMEF{15]. All these 
methods are applied to the "Parrot" image corrupted with various percentages of Salt and Pepper Noise 
(SPN). 


Table (1) presents the objective quality measurements in terms of PSNR and MSE for the "Parrot" image 
corrupted with 10%, 30%, and 50% of Salt and Pepper Noise (SPN). The results clearly indicate that the 
proposed algorithm outperforms the other filtering methods. To further analyze the effect of noise density on 
the performance, Figure (3) displays the PSNR values of the proposed algorithm and the related methods. 
These findings are visually confirmed in Figure (4), which showcases the noisy image corrupted with 40% 
SPN and the corresponding filtered images using all compared methods. Figure (4) illustrates that the 
proposed algorithm excels in noise suppression and detail preservation. 


Figure (3) clearly demonstrates that the proposed algorithm achieves superior results in both low and high 
noise densities, primarily because of the following main reasons: 


1. By employing a robust noise detection scheme 


2. The approach involves employing an adaptive window size during both the noise detection stage and 
noise filtering stage, which is determined dependent on the count of noise-free pixels within the observed 
window. 


3. Incorporating the uncertainty of fuzzy logic in both the noise detection stage and noise filtering stage. 
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Table (1): The comparative results of the proposed algorithm with different methods in terms of 
PSNR and MSE 


Noise density 
10% 30% 50% 

PSNR | MSE (x10-’) | PSNR | MSE (x10~) | PSNR | MSE (x10-’) 
Noisy image | 14.76 3.3392 10.75 8.4041 8.50 14.094 
MF [11] 28.16 0.1525 22.34 0.5829 15.07 3.1114 
HFF [12] | 27.54 0.1761 26.29 0.2346 23.90 0.4068 
FSM[13] | 35.07 0.0311 24.94 0.3201 16.54 2.2132 
AMF[I4] | 32.43 0.0571 28.78 0.1323 25.71 0.2685 
AWTMF[15] | 36.60 0.0219 31.38 0.0727 28.05 0.1564 
Proposed | 36.63 0.0217 31.74 0.0669 28.57 0.1388 


Method 


—>— Proposed 

—*— AWTMF (2013) 
—- AMF 

—e— HFF (2002) 
—— FSM (2008) 
—k— MF 


NO 


a 
22 
i 
Zz 
6 2 
a 


Noise density (SPN) 


Figure (3): Comparison chart in PSNR of the proposed algorithm with different methods 
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Figure (4): The SPN removal results of the 'Boats'' image are as follows: 
(a) Original image 
(b) Image tainted by noise contamination 40% SPN (PSNR: 9.46) 
(c) Image that has undergone filtering with the utilization of MF (PSNR: 18.34) 
(d) Image that has undergone filtering with the utilization of FSM (PSNR: 20.07) 
(e) Image that has undergone filtering with the utilization of HFF (PSNR: 25.25) 
(f) Image that has undergone filtering with the utilization of AMF (PSNR: 27.05) 
(g) Image that has undergone filtering with the utilization of AWTMEF (PSNR: 29.54) 
(h) Filtered image using the proposed algorithm (PSNR: 30.08). 


4. Conclusions 

In this study, the following key findings have been concluded: 

1. The effectiveness of fuzzy techniques in image noise reduction relies on accurately selecting the 
appropriate fuzzy sets and fuzzy rules, determining the membership function boundaries with precision, 
and employing a reliable defuzzification process. 
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The use of an adaptive window size in both the noise detection and noise filtering stages contributed to 
the superior effectiveness of the suggested algorithms in terms of objective measurements and visual 
quality, even in cases of high noise density, outperforming other methods. 


The shape of the membership functions in fuzzy set "Similar," as depicted in figures (2), is dynamically 
adjusted based on the homogeneity level of the processed window. This characteristic empowers the 
proposed algorithms to effectively differentiate between image details and noisy pixels, resulting in 
superior outcomes in relation to the elimination of noise and the preservation of image details. 
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